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Refinement of Block Motion Estimate Using Texture Mapping 



Background of the Invention 

This invention pertains to the temporal analysis of an image sequence, such as in video or 
film. In such sequences, there is substantial redundancy between successive frames. 
Most algorithms for compression of video have made good use of this temporal redim- 
dancy to reduce the data requirements for encoded video. 

In most current algorithms, images consisting of a rectangular array of pixels are divided 
into blocks, each of which is essentially a subimage consisting of a rectangle of pixels, 
often 8x8 or 16x16 pixels. These blocks may be encoded directly as in a keyframe, or 
they may make use of a prediction typically derived from temporal analysis, which is 
interframe encoding. In case of interframe encoding, the prediction makes reference to 
another frame in the sequence, often the predecessor frame, and usually provides a pixel- 
wise offset to specify a rectangular block within the prediction frame which provides a 
prediction of the particular block imder consideration in the frame which is being regen- 
erated. A residual correction to the pixel values obtained from the prediction is typically 
provided in the encoding to compensate for the inaccuracy of the prediction. 

Selection of predictive blocks for interframe coding has been an area of fertile invention 
for many years. In general, the selection of a predictive block is based on a restricted 
search of the prior frame at encoding time. Such a search typically provides an offset to a 
predictive block in whole (integral) pixel coordinates. 

Experience with video encoding and temporal analysis has revealed that viewers are very 
sensitive to small motions of objects in image sequences; viewer sensitivity is often a 
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small fraction of a pixel for large objects, such as actors in a scene. Hence many algo- 
rithms have permitted the use of predictive block offsets at subpixel resolution. Typically 
these subpixel resolutions are limited to 1/2 or 1/4 of a pixel, which allows an 
improvement in the quality of the prediction, at the cost of additional computation in 
encoding and decoding. Subpixel offsets may be obtained by interpolation of the source 
material, followed by search at higher resolution, or by analysis and interpolation of the 
error metrics which guide the search for the best integral pixel offset. 

The present invention offers a novel approach to obtaining subpixel resolution which is 
based on the use of texture mapping, a highly effective and widespread technique from 
computer graphics. In texture mapping, pixels from a reference image are mapped to a 
surface, which is arbitarily positioned and oriented in the viewing space. This permits a 
relatively low-resolution model of the geometry of the surface, complemented by a high- 
resolution representation of the visual texture of the surface. The result is efiicient and 
realistic rendering; hence the wide adoption of the technique. 

The present invention uses texture mapping to evaluate predictive offsets at arbitrary reso- 
lution. The texture mapper incorporates interpolation as necessary to provide resampling 
of the predictive image to generate the prediction. The prediction generated by the tex- 
ture mapping may be evaluated in the usual pixel-by-pixel fashion. This invention dis- 
closes a method for the use of texture mapping to successively refine a prediction offset to 
arbitary resolution. 



Summary of the Invention 



Given a soxirce image and a target image, with a block of NxM pixels in the source 
image, and a two-dimensional initial motion estimate to a matching NxM block in the tar- 
get image, the invention will refine the initial motion estimate to arbitrary precision, 
under any supplied pixel-wise metric. 

A bounding-box containing the initial motion estimate is provided, typically one pixel in 
size when the process is initiated, centered on the initial motion estimate. The error of 
the initial estimate under the supplied metric at the initial motion estimate is provided. 
The desired precision is specified as a bound on the number of refinement steps. Each 
refinement step doubles the precision of the estimate in both x and y dimensions. 

In the refinement step, the supplied boimding box is divided into four equal-sized child 
bounding boxes using the usual quad-tree subdivision, which is well known in the art. 
The center point of each child bounding box provides a new trial motion estimate. This 
new trial motion estimate is used to specify a texture map from the target image to the 
source block. The inverse image of the source block mider the texture map is typically 
not pixel-aligned to the target image. The texture map predicts the source pixels of the 
block fi-om the target image. The error of prediction is computed pixel-wise in the source 
image, using the supplied metric. If the limit on refinement steps has not been exceeded, 
the child with the smallest error is successively refined in a recursive fashion. The refine- 
ment step concludes with selection of either the initial motion estimate and its error level, 
or that of the refinement of the best of its children. 



Brief Description of Drawings 

Figure 1 is a flow diagram of the motion refinement process. 

Figure 2 depicts the source and target images as well as the source block and an initial 
pixel- wise offset to a prediction block. 

Figure 3 depicts the process of subdivision. 

Figure 4 depicts an example of three successive refinement steps. 

Figure 5 provides an altemative view of the three refinment steps in Figure 5, where the 
scale depicted is doubled with each refinement step. 

Figure 6 is a skeletal listing of the motion refinement algorithm, comparable to Figure L 

Figure 7 provides a skeletal listing of the preferred embodiment of the texture mapping 
process, using OpenGL. 
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Detailed Description of the Invention 



Figure 2 depicts the context in which the present invention operates. A source image 210 
contains a rectangular block of pixels 212. In the preferred embodiment, the rectangular 
block of pixels 212 is square, with a side length which is a power of 2 (eg, 2x2, 4x4, 8x8, 
16x16, etc). 

By any of a wide variety of block-matching means which are not disclosed here, but 
which are well known to those skilled in the art, a pixel displacement 240 has been 
obtained which relates the source pixel block 212 to a block of pixels 222 in the target 
image 220. The exact process by which such a block match is obtained is irrelevant to 
the present invention; however such a match must be provided by some means to initial- 
ize the refinement process. 

In the preferred embodiment, the source and target images 210 and 220, respectively, are 
adjacent frames in a sequence of video or film. However the process does not require 
temporal adjacency of fi-ames, and application of the invention to nonadjacent fi-ames 
may be usefiil xmder some conditions, such as where the motion between firames is very 
small. 

The source block 212 is depicted to have its lower left comer at position xq, yo in the 
source image 210. The matching target block 222 is depicted to have its lower left comer 
at postion xq + x, yo + J in the target image 220. Hence the initial displacement 240 fi-om 
source to target is given by x, y. At the outset of the refinement process, the initial dis- 
placement 240 provides the best estimate of the motion of the source block 212. 

The quality of the initial displacement as a motion estimate may be evaluated using a 
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metric which operates on pixel-by-pixel differences between the source block 212 and the 
target block 222. Various such metrics are known in the art, of which the most com- 
monly used include L\L^, and corresponding to the average of absolute pixel differ- 
ences, the square root of the average of squared pixel differences, and the maximum 
absolute pixel difference, respectively. In the preferred embodiment, the metric (aver- 
age of absolute pixel differences) is utilized, consistent with its wide use for block match- 
ing in the prior art. The selection of the error metric is independent of the present inven- 
tion, although the results of the refinement will vary as the metric is changed. 

Referring now to Figure 1, the input 100 to the process consists of an initial bounding 
box for the search 102, a best estimate of motion 104, a best error 106 associated with the 
best motion estimate 104, on which the initial bounding box 102 is centered, and a depth 
bound 108 which limits the number of refinement steps in the process. 

To start the refinement process, we obtain the initial bounding box 102 from a block 
match as depicted in Figure 2. The center of the initial bounding box is placed at x, y, 
corresponding to the pixel-wise offset determined by an external block-matching 
alorithm. In the preferred embodiment, the dimensions of the bounding box are fixed at 
1x1 pixel, hence constraining the refinement process to a bounding box ranging from 
{x-0.5,y- 0. 5) to (jc + 0. 5, >^ + 0. 5) pixels. However, the initial bounding box may be 
set to some other value, which in general will be of size Sxe. 

We obtain the initial value for the best estimate 104 from the supplied block match, 
which is presumed to provide at least a local optimum for pixel-wise matches. The best 
estimate initially is x, y, following the notation of Figure 2. 

The best error 106 is initially obtained by application of the supplied metric between the 



source block 212 and the target block 222. 



The remaining input is the depth bound 108 which limits the number of refinement steps 
and hence the resolution of the motion vector estimate. As the search space is subdivided 
by half in each direction for each refinement step, the ultimate precision of the refinement 
is related to the power of two exponent of the depth. Thus, for a depth bound of 3, the 
ultimate resolution of the refinement estimate will be 1/8 x 1/8 pixel. A depth bound of 3 
is utilized in the preferred embodiment. 

In general, the refinement will be invoked recursively, in each case provided with input 
100 comprising an initial bounding box 102 which constrains the search region for the 
refinement, a best estimate 104 and a best error 106, which may be modified if the refine- 
ment detects an improvement, and a depth bound 108 which limits the number of refine- 
ment steps and hence the precision of the resolution estimate which the refinement pro- 
cess will obtain. 

In the subdivision step 300, the initial bounding box 102 is subdivided into four child 
boundmg boxes. Turning to Figure 3, the initial bounding box 102 with center 104 at (x, 
y), occupying the region bounding by (x - SIX y - ell\ {x + 512, y + ell) is subdivided 
to four child bounding boxes, 310, 312, 314, 316, each with corresponding centers 320, 
322, 324, 326, respectively. The coordinates of the child bounding boxes and their cen- 
ters may be observed in Figure 3. 

Returning to Figure 1, in the evaluation step 400, each child boxanding box 310, 312, 314, 
316, is considered in succession. For each child bounding box, a texture mapping 450 
fi-om the target to the source is computed, utilizing the child estimate corresponding to 
the center of the child bounding box, 320, 322, 324, 326, respectively. 
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The texture mapping 450 may be performed by a variety of techniques drawn from the 
prior art. In the preferred embodiment, the freely available Mesa implemention of the 
standard computer-graphics specification OpenGL has been utilized but other implemen- 
tations could easily be substituted. Figure 7 provides a schematic of the OpenGL code 
which is used to perform the texture mapping. This is a particularly trivial usage of 
OpenGL texture mapping, utilizing a single quadrangle with two-dimensional coordi- 
nates, with a constant translation between each source, target vertex pair. 

Returning to Figure 1, the evaluation step 400 concludes with the selection of the best 
child from the foxu" children which have been evaluated. Of course, the best child is that 
which provides the smallest error imder the supplied metric. 

After evaluation, the depth boxmd is considered. In the event that the depth bound is posi- 
tive, one or more additional refinement steps remain to be performed. In this case, the 
refinement step 500 is invoked. This is a reciu*sive application of the present invention, 
however, the initial bounding box is reset to the bounding box of the best child, the initial 
error is set to the error of the best child, and the depth boxmd is decremented in this recur- 
sive invocation. Note that the best child error and corresponding estimate is reset accord- 
ing to the results of the refinement. 

Upon completion of the recursive refinement, if any, the best child error is compared to 
the best error so far. In the event that the best child error is better than the best error so 
far, the best error and the best estimate are replaced with the respective values obtained 
from the child in the reset step 600. 

In Figure 4 and Figure 5, an example of three steps of refinement is shown in two com- 
pound views. The refinement proceeds successively through SE, NW, and NE quadrants. 



corresponding to 316, 330, 340, respectively, in both Figures. In this example, the evalu- 
ation of child motion estimates at successive stages selected 326, 332, 342, respectively, 
as the best estimate from amongst the child bounding boxes. In Figure 4, successive 
views are drawn at the original scale, hence the successive search quadrants may be seen 
to shrink at each stage. In Figure 5, each successive view is drawn at double the scale, as 
indicated by the coordinates on the successive views. 

Figure 6 provides an alternative representation of the present invention, in psuedo-code. 
Due to the recursive nature of the invention, the algorithm corresponds to the refinement 
step 500. In Figure 6, the subdivision 300 from Figure 1 is carried at at 500-002. The 
loop 500-003 through 500-012 comprises the evaluation of the child estimates and the 
selection of the best child, corresponding to the evaluation step 400 in Figure L The 
optional recursive refinement step 500 is invoked at 500-014. The optional reset of the 
best estimate and error 600 from Figure 1 is carried out in 500-016 through 500-019. 



